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Abstract. We present preliminary results on the ability of self-gravitating discs to cool in response 
to their internal heating. These discs are modelled using a Smoothed Particle Hydrodynamics (SPH) 
code with radiative transfer JT|] and we investigate the ability of these discs to maintain a state of 
thermal equilibrium with their boundaries as an indication of their likelihood to fragment. 
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INTRODUCTION 

Gravitational instability has been proposed as a theory for giant planet formation |0, S], 
involving massive discs in the early stages of their evolution, and is particularly favoured 
for its rapid formation timescales. The stability of a self-gravitating disc can be described 
by the stability parameter [4], Q = where c s is the sound speed in the disc, K is the 
epicyclic frequency, which for Keplerian discs is approximately equal to the angular 
frequency, £1, E is the surface mass density and G is the gravitational constant, and 
where a value < 1 implies instability. Once the surface mass density and the rotation of 
the disc have been established, the stability is purely dependent on the disc temperature. 
Gammie [5] showed that in addition to this, a fast cooling rate is also required in order 
to overcome the stabilising compressional heating. 

The early evolution of massive self-gravitating discs has been considered using sim- 
plified cooling parameters to describe the thermodynamics [e.gjo] and using grid-based 
and hydrodynamical radiative transfer calculations [e.g. [3, |90 though results for the 
latter have differed. We use SPH with flux limited diffusion yj] to model the evolution of 
such radiative transfer discs to determine the conditions more likely for fragmentation. 



DISC SETUP 



We model a 25AU, 0.1M Q disc using 250,000 SPH particles, surrounding a 1M Q star, 
modelled using a sink particle. We use temperature and surface mass density profiles 

of T oc R 2 and E R 1 respectively, such that the aspect ratio is initially ~ 0.05. The 
temperature of the disc's vertical boundary reflects the heating from the star and this 
boundary is applied for optical depths of T < 1. It is also kept at a constant temperature 
and is initially in a state of thermal equilibrium with the midplane. 



FIGURE 1. a. Azimuthally averaged Toomre values for the initial discs (solid) and at t = 6.40RP for 
lOx (short-dashed), lx (reference; dotted), O.lx (long-dashed) & O.Olx (dot-dashed) interstellar opacities, 
b. Simulation image of a marginally stable self-gravitating disc showing high (white) to low (red) density. 

RESULTS & DISCUSSION 

The key result is that modelling discs more realistically with radiatve transfer reduces 
the likelihood of fragmentation in self-gravitating discs. Given that the boundary tem- 
perature profiles have been set, if the discs are able to respond rapidly enough to the 
internal heating, the Toomre profiles at the end of the simulations will be the same as the 
initial Toomre profiles. If the discs heat up and cannot even maintain thermal equilibrium 
with their vertical boundaries (and heat up), fragmentation is not likely. 

Figure [T^ shows that the reference and high opacity discs heat up due to internal pro- 
cesses and cannot cool fast enough to regain thermal equilibrium with their boundaries. 
However, the very low opacity disc is able to cool rapidly enough in order to maintain 
a state of thermal equilibrium with the disc boundary. This is also the case with a larger 
(300AU), cooler disc, suggesting that given more favourable (cooler) boundary condi- 
tions, these discs may potentially cool quickly enough to fragment. In fact, we find that, 
indeed, combining these properties, a larger, cooler, low opacity disc is able to fragment. 

We are currently investigating the effect of reducing the boundary temperature dur- 
ing the simulations in order to determine whether it is realistically possible to obtain 
marginally stable discs such as that shown in figure [Tt>, or to drive them to fragment. 
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